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1. Introduction 



Consider the difference equation 

C k u k -i + A k u k + B k u k+ i = Eu k , k = 1 . . . to, (f ) 

where A k , B k , C k G c mxm are complex non singular square matrices, E is a complex 
parameter, and u k G C m are unknown vectors. It is the prototype of several equations 
that occur in physics. 

At each k the equation provides u k+ \ in terms of u k and u k -i] the recursion is made 
single-term by doubling the vector and introducing the 1-step transfer matrix t k (E), of 
size 2m: 

r B k \E-A k ) 



Uk+l 




u k 







Iteration builds up the n-step transfer matrix T(E) 
vectors to steps apart: 



T(E) 



u k 

Uk-l 



t n (E) 



(2) 

■ti(E) that connects 



Ui 




U n +1 


too 







(3) 



One is often interested in the singular values oj > . . . > o 2m of T(E) (the eigenvalues 
of the positive matrix (T^T) 1 / 2 ), which describe the growth or decay of ||w n ||- The 
product of the p largest ones (p = l,...,2m) can be obtained by the formula 
oi • • -cr p = ||A p T(£)||, where (A p T)(vi A . . . A v p ) =: Tv x A 



A Tv p extends the action 



of T to antisymmetric p— forms and ||0|| is the sup norm of operators [H [2]. For real 
transfer matrices the product has the simple geometric interpretation 

Volume P{Tvx, . . . , Tv p } 



01 • • • ov 



sup 

VI... Vp 



(4) 



d2»i 



Volume P{f i, . . . , v p } 

where P{f i, . . . , v p } is the parallelogram with sides Vi G 
When the transfer matrix is the product of random matrices, Oseledets' Multiplicative 
Ergodic Theorem ensures that the singular values grow or decay exponentially in to with 
rates (Lyapunov exponents) = lim n ^oo - In a k that are independent of the realization 
[3111]. Then: 

Ai + . . . + A, 



1 



lim -ln||A p T||. 

n— >oo n 

The formula can be implemented numerically for the evaluation of Lyapunov spectra [5] . 
In the simplectic case (X m+ k = — A*) the average of the positive Lyapunov exponents is 
expressible in terms of the average distribution of eigenvalues of the Hermitian random 
matrices associated to (JTJ): 
Ai + . . . + A n 



rn 



J dE'p(E') In \E - E'\ + const. 



(5) 



The formula was obtained by Herbert, Jones and Thouless for m = 1, and by Kunz, 
Souillard, Lacroix [6] for m > 1. It is desirable to obtain similar equations for the 
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evaluation of single or other combinations of the exponents. 

In this paper the properties of a single transfer matrix are investigated. It will 
be proven that, for large n, half of its singular values have a lower bound that grows 
exponentially in n, and the other half have an upper bound that decays exponentially 
in n. Moreover, the spectrum of eigenvalues will be linked, via duality, to the spectrum 
of the difference equation ([1]) with proper boundary conditions. 

The idea of duality is simple. For a chain of length n, if Bloch boundary conditions 
(b.c.) u n+ i = e lip Ui and u = e~ ltp u n are chosen (they correspond to an infinite periodic 
chain), an eigenvalue equation is obtained: 

TIE) Ul =<?* Ul . (6) 

The condition det[T(E) — e 1Lp \ira\ = gives the nm eigenvalues E a ((p) of the difference 
equation ([I]). Then, for each eigenvalue, the whole eigenvector of the chain {u\ . . . u n ) is 
constructed by applying the 1-step transfer matrices to the initial vector (mi,m ). 
The opposite approach is also useful. The eigenvalue equation for T(E) 



Ui 




Ul 




= z 




u 




u 



T{E) 

«o 

is solved whenever (u\, . . . , w n )* is an eigenvector with eigenvalue E of the matrix 

c 2 ••. ••. 

zB„ C n A„ 



(7) 



H(z) 



(8) 



which encodes the b.c. u n+ i = zu x and u = u n /z that are implied by the eigenvalue 
equation for the transfer matrix. The statement 



Proposition 1.1 (ui, 



U r 



is a right eigenvector with eigenvalue E of the matrix 



H(z) if and only if (ui,u n / zf is a right eigenvector ofT(E) with eigenvalue z, 

translates into a determinantal identity (the duality relation, [7]) that relates the eigen- 
values of the tranfer matrix T{E) to those the associated "Hamiltonian" matrix H(z), 
that describes the difference equation of length n with generalized Bloch boundary con- 
ditions. 



It is occasionally useful to replace the parameter z with z n . The matrix H(z n ) is 
similar to the balanced matrix 

A 1 zB x |Ci 



H B (z) 



\C 2 



zB r > 



zB n _i 

ir* a 

2 w ti - rL n 



(9) 



Identities and exponential bounds 



4 




Figure 1. Left: the complex eigenvalues of a Hatano Nelson tridiagonal matrix 
(to = 1, n = 600, £ = 1) with random diagonal elements uniformly chosen in 
[—3.5, +3.5]. They lie on the line £ = X(E). The real eigenvalues correspond to states 
with localization length less than l/£. Right: the same system, with £ increasing from 
to 1 in five steps to show the expanding spectral curve. 

by the similarity relation H(z n ) = D(z)H B (z)D(z)~ 1 , where D(z) is the block diag- 
onal matrix (zl m , . . . , z n I m ). As a consequence H(z n ), H B (z) and also H B (ze lk2n ^ n ), 
k — 1 . . . n — 1, have the same eigenvalues. 

While the matrix H(z n ) remarks the value of boundary condition parameter, 

the matrix H B (z) remarks the invariance under cyclic permutations of blocks (the ring 
geometry) of the difference equation (and is numerically more tractable). 

Tridiagonal matrices of type (Q, with z = e^ real, were introduced by Hatano and 
Nelson [8] to model vortex pinning in superconductors: 

e i u k+ i + a k u k + e~ 5 w fc _i = Eu k , 

where a k are independent random entries. The model attracted a great interest as 
it gave another view of the relationship between localization and spectral response to 
b.c. variations. For zero or small £ the eigenvalues are real and all eigenvectors are 
exponentially localized with localization lengths 1/X(E). The Lyapunov exponent can 
be evaluated by Thouless' formula (JSJ), X(E) = f dEp(E) In \ E — E'\, with the average 
spectral density of the Hermitian chain (the analytic evaluation is possible in Lloyd's 
model, with Cauchy disorder [9]). By increasing £ beyond a critical value the eigenvalues 
start to gain imaginary parts and distribute along a single expanding curve [10] of 
equation £ = X(E) (see figured]). The transition has been studied also in 2D, where 
the critical value of £ for the onset of migration in the complex plane gives the inverse 
localization in the center of the band [11] . 

If the parameter £ is turned on in tridiagonal random matrices that are not Hermitian 
at the beginning, 

b k e^u k+ i + a k u k + e^c k u k -\ = Eu k , 

the phenomenon shows up differently [12]: beyond a critical value of £, an area occupied 
by the complex eigenvalues starts to be depleted, the eigenvalues being swept away and 
accumulated on an expanding "front line" of equation £ = X(E). No eigenvalues are left 
in the interior (corresponding to derealization of states) (see figure [2]). 
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Figure 2. Left: the eigenvalues of a tridiagonal matrix (m = 1, n = 800, £ = 0.5) 
with elements ak^bk,Ck chosen uniformly in [—1,1]; the "front circle" contains the 
eigenvalues that filled the circle at lower values of £. Right: the motion of eigenvalues 
(n = 100) is traced for £ changing from 0.3 to 0.6. The outer eigenvalues are numerically 
unaffected before being reached by the "front circle" . 



Though the theory presented in this paper is very general, these two models were the 
starting motivation: 

1) band random matrices have block tridiagonal structure with lower and upper 
triangular B and C matrices. Matrix elements are independent and identically 
distributed (i.i.d.) random variables. It is customary to name m as b (bandwidth 
is 2b + 1). If the probability distribution has zero mean and finite variance, and if 
n ^> b ^> 1, the spectral density of Hermitian banded matrices is Wigner's semicircle 
law, with exponentially localized eigenvectors. The localization length and its finite 
size scaling were studied numerically by Casati et al.[13j, with insight provided by the 
kicked rotor model of quantum chaos. Several properties were obtained analytically by 
supersymmetric techniques in a series of papers by Fyodorov and Mirlin |14j . 

2) Anderson model describes the propagation of a particle in a lattice with random site 
potential. After choosing a (long) direction of length n, the diagonal blocks = T + D^ 
describe the sections of the lattice with m sites each (T is the Laplacian matrix for the 
transverse slice and is a random diagonal matrix with i.i.d. elements). The hopping 
among neighboring slices is fixed by Bk = Ck = I m . The random site-potential is 
usually chosen uniformly distributed in [— w/2,w/2] {w is the disorder parameter) (the 
literature is vast, see [15] for a mathematical introduction). 

In both models the transfer matrix is a product of random matrices and, for n — > oo, 
it provides a non random Lyapunov spectrum [T6| [T7J [18] . The inverse of the smallest 
Lyapunov exponent is the localization length. 

Localization affects the response of energy values to variations of b.c. [HI [20] . This dual 
way of viewing localization: through decay of eigenvectors (transfer matrix) or response 
of energy levels to b.c. variations (Hamiltonian matrix), is hidden in the duality identity 
among the eigenvalues of T(E) and of H(z). 
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The first two sections provide algebraic properties that relate a generic transfer 
matrix to its Hamiltonian matrix. Some of them appeared in previous papers, but 
receive here a consistent presentation. In particular, they are the spectral duality and 
the expression of T(E) in terms of the resolvent of the Hamiltonian matrix with open 
b.c. 

Next, a theorem by Demko, Moss and Smith [2TJ on the decay of matrix elements of the 
inverse of a banded matrix is presented. It is used here to prove that a 2m x 2m transfer 
matrix has m singular values growing exponentially with the length of the chain, and 
m singular values decaying exponentially. This new result reflects on a single matrix a 
property of random matrix products. 

The rest of the paper deals with identities; duality and Jensen's identity give an 
expression for the exponents £ a = i In \z a \, where z a are the eigenvalues of the transfer 
matrix, in terms of the eigenvalues of the associated matrix H(z). Hadamard's inequality 
for determinants of positive matrices supports the idea that the eigenvalues z a have a 
leading exponential growth in n. The discussion of the relevant case of Hermitian 
difference equation ends the paper. 



2. Transfer matrix and duality 



Some general facts about transfer matrices are presented. By construction T(E) is a 
polynomial in E of degree n, T(E) = E n T n + . . . + ET\ + T , with matrix coefficients. 
However, its determinant is independent of E: 

This implies that T{E)~ X is again a matrix polynomial in E [22]. Actually T(E)~ l is 
similar to the tranfer matrix of the inverted chain. Let's introduce the two matrices of 
inversion, of size 2m x 2m and nm x nm: 



o-. 



l m 

l m o 



J 



Ir. 



Proposition 2.1 Let T(E) be a transfer matrix and H(z) the associated matrix, and 
let T(E) J be the transfer matrix associated to H J (z) = JH(z)J (the inverted chain); 
then: T(E)- 1 = a x T(E) J a x . 

Proof: T(E)- 1 = [t n (E) ■ --t^E)}- 1 = t^E)' 1 ■ ■ -t n (E)~ l . The combination 

C k \E-A k ) -C?B M 

lm. 



gives the structure of a 1-step transfer matrix. Multiplication yields the result. 



□ 
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Proposition 2.2 In the expansion of the characteristic polynomial of the transfer 
matrix, 

det [zl 2m - T{E)\ = z 2m + . . . + a k (E)z 2m - k + ... + a 2m „ k (E)z k + ... + a 2m , 

the coefficients a k (E) and a 2m _ k (E) are (in general different) polynomials in E of degree 
kn (k = 0, m ). 

Proof: Let z±, . . . , z 2m be the eigenvalues of T(E). The coefficients 
a k = (-l) k ^ z h ---z ik , k = l...m, 

i\<...<i k 

can be expressed as combination of traces of powers of T(E) of degree k: a\ = — tr T(E), 
a 2 = |[tr T{E)} 2 — | tr [T(E)% etc. Since T(E) = E n T n + ...+T , the coefficient a k is a 
polynomial of degree kn in E. The remaining coefficients a 2m - k are discussed differently. 
The point is that a 2m — Z\ • • ■ z 2m = detT(E) is independent of E and the coefficients 
can be written as 

a 2m -k = (~l) fc ^2 Z ^ ' ' ' Zi 2m-k — (— l) fc «2m 22 ( Zi i'" Z ik) 1 

<...<it 

Therefore, a 2m -i = -a2m tr[T(E)~ 1 ], a 2m -2 = a 2m |[tr T(E) -1 ] 2 - a 2m | tr[T(E)~ 2 ], etc. 
Since also T(E)~ 1 is a polynomial matrix of degree n in i£, a 2m _ k is a polynomial of 
degree kn in E. □ 

Theorem 2.3 (Duality) 

^-^H-.r ^ (id 

Proof: According to proposition 12.21 the leading term in the expansion in i£ of 
det [0l 2 m — T(E)} coincides with the leading term in the expansion of det[2;I 2m — E n T n ], 
which is (-z) m E nm det[ J B 1 --- J B n ]- 1 . The leading term of det[M 2m - H{z)\ is E nm . 
Since by proposition 11.11 the two polynomials, for given z, have the same zeros in E, 
they must be proportional by a constant. □ 

This relation among characteristic polynomials is a "duality identity" as it 
exchanges the roles of the parameters z and E among the two matrices: z is an eigenvalue 
of T(E) if and only if E is an eigenvalue of the block tridiagonal matrix H(z). I gave 
different proofs of it [TJ |23l El]. With z — 1 it is a tool for computing determinants of 
block tridiagonal or banded matrices with corners. 

The eigenvalues of H{z) make the l.h.s. of duality equal to zero, i.e. there is at least a 
complex factor Zi(E) — z = 0. This means that an eigenvalue E is at the intersection of 
a line (^(-E 1 )] = \z\ and arg Z{(E) = arg z. By changing only the parameter arg z, the 
eigenvalues move along spectral lines \zi(E)\ = \z\. For tridiagonal matrices (m = 1) 
there is a single spectral curve (figured]), for m > 1 several spectral curves appear [25] 
(see figure [3]) . 
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Figure 3. The eigenvalues of the Hamiltonian matrix (n = 8, m = 3) of a 2D Anderson 
model on a lattice 3x8, diagonal disorder parameter w — 7. The b.c. parameter is 
z = cxp(n£ + i(p) with £ = 1.5. As (p changes, the 24 eigenvalues (of the balanced 
matrix) trace m = 3 closed loops of equation In |zfe(i?)| = 1.5. 



A more symmetric duality relation results from multiplication of the dual identities 
for (T — z) and (T — 1/z): 



det 



T(E) + T(E)- 1 - [z+- z ] Io„ ; 



det[M wm - H{z)\ det[M wm - H{l/z)} 
det[B l ---B n ] dettd--.Cn] 



3. Transfer matrix and resolvent 



Equation (CQ) with open b.c. uo 
matrix 



and u n+ i = 0, is the eigenvalue equation for the 



h 



C 2 



Cr, 



y4„ 



(12) 



Let (ui, . . . ,u n Y be a (right) eigenvector of h with eigenvalue £7; then u\ and u n are 
both nonzero, or the whole vector would be null by the chain recursion. With the block 
partition 



T(E) 



T(E) 1A T(E) li2 
T{E) 2>X T(E) 2>2 



03]) gives T(E)i t iUi = and T(E) 2 ^ui = u n . This means that detT(£')i j i = whenever 
det[-EI nm — h] — (and detT(i?) 2j i ^ 0). The following duality relation results: 

Proposition 3.1 (duality for the open chain) 

det[M nm - h) = detT(E) ljl det^ - - ■ B n ] (13) 

Proof: by construction T(E)u = E n (B\ ■ ■ ■ B n )~ 1 + lower powers in E. Then both 
det[EI nm — h] and det T(E)n are polynomials in E of degree nm. Having the same 
roots, they are proportional. □ 
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The blocks T(E)i 2 and T(E) 21 are polynomial matrices of degree n — 1 in E, and 
T(E) 22 has degree n — 2. The four blocks can be evaluated in terms of the corner blocks 
of the resolvent matrix 



g(E) = [ft - EI nn 



9l,r 



9n,l ' ' ' Qn^n 

The corner matrices C\ and B n are absent in ft but enter in the definition of T(E) 
through the 1-step factors t\(E) and t n (E), and will be accounted for. 

Proposition 3.2 Let g Lj G C mxm (a,b = 1 . . .n) be the blocks ofg(E). Then 



T(E) 



-B-\g 

9n,n(9l,r 



l,n) 



i-l 



9r. 



-B- 1 (g ljn )- 1 g ljl C 1 

) fi'l.lCl ~~ 9n,\C\ 



(14) 



Proof: write the identity [ft, — £'I nm ]^(£') = I nm for the block indices i — 2 . . . n — 1 and 
= l,n: Cig^i t k + (Ai — EI m )gik + Big i+ i t k = 0. The recursive relations are solved by 
the transfer matrix method and give a matrix relation among the corner blocks: 



9n,l 9n 
9n~l,l 9n-l,r. 



t n -i{E) ■ ■ -t 2 (E) 



92,1 92,r. 
91,1 9l,r. 



Left multiply both sides by t n (E) and simplify l.h.s. by means of the identity 
C n g n -x,k + {A n - EI m )g nik = 4,„I m . Insert ti^t^Ey 1 = l 2m in the r.h.s. to obtain 
T(E)t~l l , and simplify the action of if 1 by means of the identity {A\— EI m )gi jk +B 1 g 2jk = 
Si fcl m . The useful factorization is obtained: 





9n,l 



9n,n 



-1 



T{E) 



91,1 9l,n 

-cy 



(15) 
□ 



A matrix inversion and multiplication give the result. 

The representation provides the transfer matrix through a large matrix inversion, rather 
than multiplications. It was employed by Kramer and MacKinnon [26J in a numerical 
proof of one-parameter scaling for the localization length of Anderson's model. 



4. Exponential inequalities 

Products of random matrices are known to exhibit Lyapunov exponents that are 
asymptotically stable and self-averaging, i.e. independent of the length n and of the 
realization of the random product. In the present deterministic approach a single chain 
is considered, and it will be shown that it is possible to give exponential bounds on the 
eigenvalues for long chains, that justify the introduction of exponents. 

Demko, Moss and Smith [21] made the very general statement that, loosely speak- 
ing, the matrix elements of the inverse of block tridiagonal or banded matrices decay 
exponentially from the diagonal (see also [271 EH])- I here present their interesting proof 
adapted to the block partitioning of matrices. I then apply it to the matrix g(E) to 
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obtain bounds for the singular values of T(E). 



The main ingredient is the best approximation of the function (x — a)" 1 on the 
interval [—1, 1] (|a| > 1) by a polynomial of degree k, which was obtained by Chebyshev 
together with the determination of the error [29]. With proper rescaling it is |21j : 



Lemma 4.1 Let Pk be the set of real monic polynomials of degree k, let [a,b] be an 
interval of the positive real line, with a > 0. Then: 



inf 

PfcGPfe 



sup 

x£[a,b] 



Pk{x) 



X 



Cq 



k+l 



c 



(Vb + y/a)' 



Vb 



(16) 



(17) 



2a6 y/b + yfa 

If A is a block tridiagonal matrix with blocks of size m x m and if Pk{x) is a polynomial 
of degree k, the blocks pk{A) it j of the matrix Pk{A) are null for \i — j\ > k. 
Let A be a positive definite block tridiagonal matrix, with inverse A' 1 . If j] 
denotes any matrix element in the block (A" 1 )^ then, for any monic real polynomial of 
degree k = \i — j\ — 1, it is: 

\A-%j]\ = \A-%j}-p k (A)[z,j]\ 

< \\A- 1 -pk{A)\\ = sup ' 



xesp(A) 



A 



< sup 



where \\A\\ = supu x u =l \\Ax\\ is the operator nor and the spectral theorem is used. 
In the last line [a, b] is the smallest interval containing the spectrum of eigenvalues 
sp(A). Next, the inf is taken over the polynomials pk- The lemma states that the 
minimum exists, and the error gives the main inequality. Note that for \i — j\ = 0: 
I^MI < WA^W = l/a. Therefore: 

Theorem 4.2 (Demko, Moss and Smith) Let A be a positive definite block 
tridiagonal matrix, with square blocks of size m, let [a, b] be the smallest interval 
containing the spectrum of A, let v4 -1 [z, j] be any matrix element in the block (A^ 1 )^. 
Then: 

Cq li ~ jl for |i — jl > 1 

(18) 

l/a for i = j 

where q < 1 and C are specified by eg. (fT7|). 

Demko et al. also proved an extension of the theorem to a matrix A that is block 
tridiagonal invertible but fails to be positive. An estimate for A" 1 is obtained by noting 
that A~ l = A> (AA')~ X . The matrix AJV is block 5-diagonal positive definite, and a 



\A-%3]\ < 



| For any matrix A with matrix elements A rs it is \A r 
canonical unit vectors and Schwarz's inequality is used. 



(e r |^4e s )| < ||-Ae s || < \\A\\, where e» are 
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polynomial pk(AA^) is a matrix whose blocks are null if \i — j\ > 2k. The previous 
theorem applies, with [a, b] being the smallest positive interval containing sp(AA'): 

\(AAl)-%j]\<-j=q&-*, \i-j\>2 

The extension of the theorem is here written in the block notation, with minor changes 
from the original paper: 

Theorem 4.3 Let A be an invertible block tridiagonal matrix with square blocks of size 
m, let [a,b] be the smallest interval containing sp(A^A), let be any matrix 

element in the block (A^ 1 )^. Then: 

\A-%j}\<C iq ^-i\ (19) 
d = - (||A- M || + II A,, || + HA+ull) , (20) 

Q 

where q < 1 and C are given in ( (77| ). 
Proof: in terms of block multiplication: 

(A-% = (i4 t ) i , i _ 1 [(AA + )~ 1 ]i_u + (A t ) l> -[(AA t )~ 1 ] i)i + (A ] ) iMl [(AA ] y l ] i+ld . 

The sup norm, the triangle inequality, the property ||AB|| < ||A||||5||, and the bound 
on (AA) _1 give: 

\\(A-%\\ < HA-LilUKAAt)- 1 ^!! + WA^WUAA^-^J + \\A l+1 4\\(AA^-\ +1 J 

< (llA-dk^- 1 ! + ||A M ||g^^i + ||A+ M ||?5i^' +1 i) 

< 5(1)^11 + 11^11 + 11^-11) 

q 

If j] is any matrix element in the block (A^ 1 )^, it is |A _1 [i, j]\ < ||Aj||- D 

Given an invertible matrix A, the condition number of A is [2J: 

cond (A) =: ||A|| || ^s.— 1 1| 

In general it is cond (A) > 1. If a and 6 are the extrema of the spectrum of a positive 
matrix P it is b = \\P\\ and 1/a = ||P _1 ||; then b/a = cond (P). 

Since ||A(4J|| = ||^4|| 2 , it is cond (AA^) = [cond(A)] 2 and the parameters in theorem 14.31 
are: 

_ cond(A) - 1 _ (cond(A) + l) 2 

g ~cond(A) + l' 2||A||2 



Theorem 14.31 is applied to the corner blocks of the resolvent g(E) — [h — EI nm \~ x , 
E ^ sp(h), which enter in the representation (j!4p of the transfer matrix. The number 
cond (h — E) defines the parameters q < 1 and C. 
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Proposition 4.4 If g[l,n] and g[n, 1] are matrix elements of the corner blocks g\ n and 
9ni of g(Et), then the following inequalities hold: 

\g[l,n}\ < Cdl^-^l + II^IDg^- 3 ), (22) 
|<7[n,l]| < C(\\A n -E\\ + \\C n \\)q^ n -V (23) 

where Ai, B\, A n , B n are the blocks in the first and last row of h. 

We prepare for the main theorem with the following lemma: 

Lemma 4.5 The singular values 9 k of the block T u ofT(E) are exponentially large in 
n: 6 k > q- n / 2 /K. 

Proof: from (Tu)" 1 = -g ln B n it follows that: tr[(T 1 t 1 T 11 )^ 1 ] = tr[fl n £Mn0irJ < 

2 "i? nJ Bt||||^ ln || =m 2|| Sn ||2||^ n ||2< m 2|| Bn ||2 C 2 ( || Ai + 11^11)2^-3 = . R 2 ? n 



m 

Since tr[(T 1 \T 11 ) _1 ] = Y^k=i @k 2 ■> ^ turns out that each singular value of T n is larger 
than q- n / 2 /K. □ 

Main Theorem 4.1 If q < 1 and n is large, the transfer matrix T(E) has m singular 
values larger than -j^q~ n ^ 2 and m singular values smaller than Kq n l 2 . 
Proof: Let 9i > . . . > 8 m be the singular values of the block T n , and let a± > . . . > cr 2m 
be the singular values of T(E). The interlacing property (Theorem 7.12 of ref. |30j) 
states that: 

Vk>6k> o- m+ k, k = l,...,m 

Therefore, there are at least m singular values of T(E) that are larger than -j^q~ n ^ 2 - 
Since the same conclusion holds true for T(E)^ 1 , which is similar to a transfer matrix 
by proposition I2.1[ there are precisely m singular values of T(E) that are larger than 
j?q~ n ' 2 , and m that are smaller than Kq n l 2 . □ 

5. Jensen's formula and the exponents 

The two sides of the duality relation are determinantal expressions of the same 
polynomial in two variables, F(z, E) =: det [zl2 m ~ T{E)]. Let z±, . . . , z 2rn be the zeros 
in the variable z (the eigenvalues of T(E)) with \zi\ > ... > \z2 m \, and let Ei, . . . , E nm 
be the zeros in E (the eigenvalues of H(z)). It is convenient to introduce the exponents 
of the transfer matrix: 

£ fc =: - In \z k \ 
n 

Remark 5.1 The exponents are not to be confused with the Lyapunov exponents, which 
are defined in terms of the positive eigenvalues x k of the matrix T^T . It has been shown 
( with less general T ) that also T^T is the transfer matrix of a block tridiagonal matrix, 
so the same discussion may be applied to them 



The sum of the exponents is - In | det T(E)\. Then: 

2m . n 

& = - Yl ( ln I det C i\ ~ ln I det B i\) ( 24 ) 



n 

k=l j=l 
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Some general analytic results are now given, based on the following theorem of complex 
anaysis [31 J: 



Theorem 5.2 (Jensen) If f is holomorphic and /(0) 7^ ; and z\ . . . z n are its zeros 
in the disk of radius r, then: ^ In |/(re ie )| = In |/(0)| — J2 k ^ n (\ z k\/ r )- 

The theorem is applied to F(z, E) as a function of z, resulting in a relation between a 
sum of the exponents and the spectrum of the Hamiltonian matrix [32J: 

Proposition 5.3 

= ^Tnt% ln l det ^( ex P[< + M) - E]\ - ^ £ ln l det ^1 ( 25 ) 

Proof: Jensen's theorem with z = e n ^ +ld gives in the r.h.s. the sum of exponents 
contained the disk of radius e n ^\ 

/ - In I F(e<+ ie , E) \ = ln | det T(E) | + n - t- h )0((- - 
Jo Z7T k=i 

The dual expression is used in the l.h.s.: ln \F\ = mn£ + ln | det[H (e n ^ +1<p ) — E]\ — 
^■InldetSj-l. □ 

A derivative in the variable £ of ( 1231) gives the counting functions of exponents 
N(£, E) = — £,a(E)), which is also obtainable by Euler's formula for the zeros Zk 

of the F(z,E) [35]. 

Hadamard-Fisher's inequality [21 [30] states that if Mi,...,M n are the diagonal 
blocks of the positive matrix A* A, then | det A\ 2 < det Mi ■ ■ • det M n . 
The inequality is applied to the r.h.s. in eq. (l2"5"j) . with the balanced matrix H B (e^ +1 ^ n ): 

2m ^ n 

]T(£ - - &) - m£ < — J> |det (26) 

fe=i n j=i 

1 n 

+^ E ln det [(4 - ^) - E ) + ^ B l B * + e-*cjc* 

fe=l 

If the norms of matrices Ai Bi and Cj are bounded by some constant for all i, and m is 
fixed, the sum in l.h.s. of inequality remains finite for any length n, as the r.h.s. is an 
average value for the blocks. 



Corollary 5.4 The sum of the positive exponents is obtained from /[2o\) with £ = and 
by means of eq. [Efy 

2m „27r in i 

V6%) = - / — ln|det[i/(e ie )- J E;]|--ln|det[ J B 1 --- J B n ]| (27) 
k=1 n J An n 
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The identity is exact and applies to a single transfer matrix. It is reminiscent of the 
formula (jSJ) for the sum of the Lyapunov exponents of random transfer matrices. The 
"angular average" replaces the ensemble averaged density of eigenvalues p(E), which 
was extended to tridiagonal non-Hermitian matrices in [33| 134] . 



6. The Hermitian difference equation 



Most of the literature concentrates on the Hermitian case. However, as duality requires 
z to be a complex parameter, the matrix H(z) fails to be Hermitian unless \z\ = 1; 

A 1 B x \Bl 

B\ ■: 

B n -i 



H(z) 



zB n 



B 



n-l 



A k = At 



(28) 



A useful symplectic property holds for the transfer matrix (in transport problems it 
describes flux conservation, [7]), and implies that exponents come in pairs ±£ a : 

-Bl 



Proposition 6.1 



T(E)^ n T(E) = S r 





B n 

Proof: in the factorization T(E) = t n (E) ■ ■ -ti(E), the factors th{E) (k = 2 . . . n) have 
the property tk{E)^T,ktk{E) = E^-i- The factor t\ that contains the boundary blocks, 
closes the loop: ti(EY'Eiti(E) = S n . □ 

Corollary 6.2 If E is real, the eigenvalues ofT(E) different from ±1 come in pairs z, 
1/z. The associated exponents are opposite. 

Proof: If T(E)u = zu, the symplectic property implies that T(EyT: n u = l/zT, n u i.e. 
1/z is an eigenvalue of T(E). Moreover, if \z\ ^ 1, then u^S n n = 0. 

Proposition 6.3 If Im E ^ then T(E) has no eigenvalues on the unit circle. 
Proof: for Im E ^ and z = e ld it is always det[i? — H(e ie )} ^ because H(e ie ) is 
Hermitian and it has real eigenvalues. Therefore, by duality, det[T(£') — e l6, I 2 m] never 
vanishes. □ 

A degeneracy occurs in the exponents of the real transfer matrix of a real symmetric 
difference equation (the Anderson model is a notable example, but remind remark I5TTT) : 

Proposition 6.4 Let the matrices Ak be real symmetric and Bk be real invertible. For 
SeK, the real eigenvalues ofT(E) come in pairs z, 1/z, the complex ones also have 
the conjugated pair 1, 1 /z. 

Proof: if z is a complex eigenvalue of T(E) not in the unit circle, then also ~z, l/~z and 
1/z are distinct eigenvalues, and exponents are doubly degenerate opposite pairs. If z 
is a real eigenvalue, then 1/z is an eigenvalue. Therefore an eigenvalue (real or not) is 
always paired to the eigenvalue 1/z. □ 
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